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The cross sections for D D and D D~ production at 3.773 GeV have been measured with BES- 
II detector at BEPC. These measurements are made by analyzing a data sample of about 17.3 pb -1 
collected at the center-of-mass energy of 3.773 GeV. Observed cross sections for the charm pair 
production are radiatively corrected to obtain the tree level cross section for DD production. A 
measurement of the total tree level hadronic cross section is obtained from the tree level DD cross 
section and an extrapolation of the R u d s below the open charm threshold. 

PACS numbers: 



I. INTRODUCTION 

Around the center-of-mass energy of 3.770 GeV, the 



open charm pairs, D°D° and D + D~, are mainly pro- 
duced in ?/;(3770) decays. So the measurement of the 
cross sections for D°D° and D + D~ production at an en- 
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ergy point around 3.770 GeV is very important in the 
understanding of V(3770) decays. A coupled-channel 
model 0] predicts that the cross section for ^(3770) pro- 
duction is about 3 nb and that the ?/>(3770) decay exclu- 
sively into D°D° and D + D~ . Experimental results on 
the measurement of the D°D°, D + D~ and DD cross sec- 
tions can be used to test the theoretical prediction. Also 
the measured values of the cross sections can be used to 
determine the branching fraction for ^(3770) — > nonDD 
using the measured cross section for ^(3770) — * hadrons 
at the same energy point. The determination of the par- 
tial width of ■0(3770) — * nonDD has great interest since 
it would be helpful for investigating the mixing between 
S and D waves in its wave function 0,and in turn to 
help in developing the Potential Model [l| . 

In addition, by adding the tree level open charm cross 
section to an extrapolation of the tree level hadronic 
cross section for the light hadron production in the re- 
gion below the open charm threshold, the total tree level 
hadronic cross section can be obtained [30- The tree 
level cross section for inclusive hadronic event produc- 
tion in the e + e~ annihilation at all energies is needed 
to calculate the effects of vacuum polarization on the 
parameters of the Standard Model. The largest uncer- 
tainty in this calculation arises from the uncertainties 
in the measured inclusive hadronic cross sections in the 
open charm threshold region. Traditionally, the tree level 
hadronic cross sections are measured by counting inclu- 
sive hadronic events. Using the measured cross sections 
for the DD production in the charm threshold region, we 
can also measure the tree level inclusive hadronic cross 
sections. 



II. THE BES-II DETECTOR 

BES-II is a conventional cylindrical magnetic detector 
that is described in detail in Ref. [5j. A 12-layer ver- 
tex chamber (VC) surrounding the beryllium beam pipe 
provides input to the event trigger, as well as coordinate 
information. A forty-layer main drift chamber (MDC) 
located just outside the VC yields precise measurements 
of charged particle trajectories with a solid angle cov- 
erage of 85% of 47r; it also provides ionization energy 
loss (dE/dx) measurements which are used for particle 
identification. Momentum resolution of 1.7%yl + p 2 (p 
in GeV/c) and dE/dx resolution of 8.5% for Bhabha 
scattering electrons are obtained for the data taken at 
y/s = 3.773 GeV. An array of 48 scintillation counters 
surrounding the MDC measures the time of flight (TOF) 
of charged particles with a resolution of about 180 ps for 
electrons. Outside the TOF, a 12 radiation length, lead- 
gas barrel shower counter (BSC), operating in limited 
streamer mode, measures the energies of electrons and 
photons over 80% of the total solid angle with an energy 
resolution of cr E /E = 0.22/VE (E in GeV) and spatial 
resolutions of er^ = 7.9 mrad and oz = 2.3 cm for elec- 
trons. A solenoidal ma,enet outside the BSC nrovides a 



0.4 T magnetic field in the central tracking region of the 
detector. Three double-layer muon counters instrument 
the magnet flux return and serve to identify muons with 
momentum greater than 500 MeV/c. They cover 68% of 
the total solid angle. 



III. DATA SAMPLE AND METHOD TO 
DETERMINE a DD 

The data used for this analysis were collected at the 
center-of-mass energy of 3.773 GeV with the Beijing 
Spectrometer [5j at Beijing Electron Positron Collider. 
The total integrated luminosity of the data set is 17.3 
pb _1 , which is obtained based on analysis of large angle 
Bhabha scattering from the same data set. 

The measurements of the cross sections for the D°D° 
and D + D~ production are made based on the analysis of 
singly tagged D° and D + events. At the center-of-mass 
energy y/s = 3.773 GeV, the D° (through this paper, 
charge conjugation is implied) and D + are produced in 
pair via the process of 



e + e 



D°D°,D^D' 



(1) 



The total observed number N D o (N n + ) of D° (D+) 

ta S tag 

meson and the observed cross section er^ s 50 (cr^ D _) 
are related as 



and 



obs_ 
>D°D° 



obs 
T D+D- 



tag 



2 x L x B x e' 



N 



d: 



2 x L x B x e 



(2) 



(3) 



where L is the integrated luminosity of the data set used 
in the analysis, B is the branching fraction for decay 
mode in question, and e is the efficiency determined from 
Monte Carlo for reconstruction of this decay mode. 

In the measurements of the cross sections, the singly 
tagged neutral and charged D mesons are observed in the 
invariant mass spectra of the daughter particles from the 
D decay. 



IV. DATA ANALYSIS 



A. Event selection 



The neutral and charged D mesons are reconstructed 
in the final states of K ~tt + , K~Tr + ir + ir~ and K~ir + ir + . 
Events which contain at least two reconstructed charged- 
tracks with good helix fits are selected. In order to ensure 
well-measured 3-momentum vectors and reliable charged- 
particle identification, the charged tracks used in the sin- 
gle tag analysis are required to be within \cos9\ <0.85, 
where is the axisrle with resnect to beam direction. All 
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TABLE I: Singly tagged D° and D+ samples. Where the 
Mat is the fitted mass of singly tagged D meson, the N?% s 
is the observed number of singly tagged D meson and the 
ND tag is the "true" number of the singly tagged D meson 
after correcting the contamination from other decay modes. 



Tag Mode 


Mm [MeV/c 2 ] 


TV" 1 - 

"tag 


N Dtas 




1865.5 ±0.1 


1642.8 ± 49.9 


1627 A ± 49.9 




1865.4 ±0.1 


1327.2 ±48.6 


1299.1 ±48.6 


K--K+-R+ 


1870.2 ± 0.1 


2029.3 ± 57.4 


2010.8 ± 57.4 



tracks must originate from the interaction region, which 
require that the closest approach of the charged track in 
xy plane is less than 2.0 cm and the z position of the 
charged track is less than 20.0 cm. Pions and kaons are 
identified by means of TOF and dE/dx measurements. 
The pion identification requires a consistency with the 
pion hypothesis at a confidence level greater than 0.1%. 
In order to reduce misidentification, the kaon candidate 
is required to have a larger confidence level for a kaon 
hypothesis than that for a pion hypothesis. 



B. Analysis of inclusive D meson events 




Fitted Mass (KV) (GeV) 




t I I I I I r"— i- 1 

1.82 1.83 1.84 1.85 1.86 1.87 1.88 1.89 

Fitted Mass (K'jiYjT) (GeV) 




1.82 1.83 1.84 1.85 1.86 1.87 1.88 1.89 



Fitted Mass (K'jtY) (GeV) 



Taking advantage of the DD pair production, we use 
a kinematic fit to candidate D° or D + decays to improve 
the ratio of signal to noise and mass resolution in the in- 
variant mass spectrum. The energy-constraint is imposed 
on the measured momenta of the D daughter particles via 
the kinematic fit to improve the measured charged track 
three-momenta. Events with a kinematic fit probability 
greater than 1% are accepted. If more than one combi- 
nation satisfies the fit probability greater than 1%, the 
combination with the largest fit probability is retained. 

The resulting distributions of the fitted masses of Kmn 
(to = 1, or 2, or 3) combinations, which are calculated 
using the fitted momentum vectors from the kinematic 
fit, are shown in figure 1. The signals for D° and D + 
production are clearly observed in the fitted mass spec- 
tra. A maximum likelihood fit is performed to the mass 
spectrum, a Gaussian is used to model the signal shape 
and a special function [6] is used to describe the back- 
ground shape. The D° and D + yields obtained from this 
fit is given in table I. 

The same Kmix combinations from other decay modes 
can also pass the above selection criteria and the distri- 
butions of the fitted masses of the combinations are with 
a small peak around the masses of D mesons. The rates 
of the contaminations are evaluated to be 0.0094, 0.0212 
and 0.0091 for the D° -> R-TT+, D° -> R-n+n+Tr- and 
D + — > K~n + ir + , respectively. After correcting the ob- 
served numbers of the singly tagged D events for these 
combinations, the "true" numbers of the D signal events 
for the three sinelv taffsred D modes are obtained to be 



FIG. 1: Distribution of the fitted masses of the Kmir (m=l, 
or 2, or 3) combinations for three singly tagged modes, where 
figure (a) and figure (b) are for the decay modes of D° — > 
K~ty + and D° — > K~tv + tt + tt~ , respectively, and the figure 
(c) is for the decay mode of D + — » K~it + -k + . 



1627 ± 50, 1299 ± 49 and 2011 ± 57. Table I summarizes 
the results of the inclusive D analysis. 



V. OBSERVED CROSS SECTIONS FOR D°D° 
AND D + D~ PRODUCTION 

A. Monte Carlo Efficiency 

To estimate the reconstructed efficiencies of D° — > 
K-JT+, D° -> A^tt+tt+tt- and D+ ->■ K-tt+tt+, the 
Monte Carlo samples of DD production and decays are 
generated according to equation (1), where the ratio 
of the neutral over the total DD production cross sec- 
tion is set to be 0.58. Both D and D mesons decay 
to all possible modes according to the branching frac- 
tions quoted from PDG 7J. The generated events are 
simulated with a GEANT based Monte Carlo package. 
All decay processes which contribute to the decay modes 
in question are considered in estimating the efficiencies. 
Detailed Monte Carlo studies give the efficiencies to be 
(35.26 ±0.19)%, (13.73 ±0.09)% and (25.00 ±0.13)% for 
the reconstruction of D° — ► K~tt + . D° — ► A'~7r + 7r + 7T - 
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TABLE II: Summary of the observed cross section times 
branching fraction. 



TABLE IV: Comparison of the observed cross section with 
that measured by MARK-III [Toll experiment. 



Mode 



a DD x B [rib] 



D 







D° -> K-tt+tt+tt- 
D+ -> K--K+-K+ 



0.133 ±0.004 ±0.008 
0.273 ±0.010 ±0.025 
0.233 ± 0.007 ±0.018 



<s [nb] <rSS [»b] 

(This experiment) (MARK-III) 

g em = 3.773 GeV £ cm = 3.768 GeV 

cr D o B o 3.58 ±0.09 ±0.31 2.90 ±0.25 ±0.30 
a D+D - 2.56 ±0.08 ±0.26 2.10 ± 0.30 ± 0.15 



TABLE III: A comparison of <jd x B measured by this exper- 
iment, MARK-I and MARK-II experiments. 



a D X B [nb] a D X B [nb] a D X B [nb] 

Mode (This experiment) (MARK-II) (MARK-I) 

Ecm = 3.773 GeV E cm = 3.771 GeV E cm = 3.774 GeV 

K~n + 0.27 ±0.02 0.24 ± 0.02 0.25 ± 0.05 

K-tt+tt+tt~ 0.55 ± 0.05 0.68 ±0.11 0.36 ± 0.10 

K~ 7T+7T+ 0.47 ±0.04 0.38 ±0.05 0.36 ± 0.06 



and D + — > if 7r + 7r + decay modes, respectively. 

B. Observed Cross Sections 

Inserting the number of the singly tagged D events 
and the efficiencies for each of the three decay modes, 
the <t d] j x B of D°D° and D + D~ are obtained and the 
results are shown in table II. The first error is statis- 
tical and second systematic which arise from the un- 
certainty in the measured luminosity (~ 3%), tracking 
(~ 2% per track), particle identification (~ 0.5%/track), 
kinematic fit (~ 1%), fitting parameters (~ 3%) and 
Monte Carlo statistics (~ 0.6%). The total systematic 
uncertainty is obtained by adding all systematic uncer- 
tainties in quadrature. A comparison of our measured 
<7nx B with that measured by MARK-II and MARK- 
I is given in Table III. The observed cross section for 
D + D~ production is obtained by dividing the <j d q x B 
by branching fraction quoted from PDG0, which gives, 

a°$. D - = (2.56 ± 0.08 ± 0.26) nb. (4) 

The observed cross sections for D°D° production are 
a°^- = (3.50 ±0.11) nb and o°$ & = (3.66 ±0.13) nb, 
which are determined from the analysis of the singly 
tagged modes of D° — > K~ir + and D° — > K~ir + ir + ir~ 
respectively; where the errors are statistical. Averag- 
ing the two observed cross sections for D°D° production 
gives the average of the observed cross section for D°D° 
production to be 

a Doi5o = (3-58 ± 0.09 ± 0.31) nb, (5) 

where the first error is statistical and second systematic 
which is estimated based on the avcraecd three charged 



tracks in the two modes of the single tags. Adding the 
observed cross sections of the neutral and charged modes 
together gives the observed cross section for DD produc- 
tion to be 

o-° b £ = (6.14 ±0.12 ±0.50) nb, (6) 

where the first error is statistical and the second system- 
atic. In the estimation of the systematic error of the 
p^fi, sources of systematic uncertainty are segregated 
into components that are common or independent for D 
and D + measurements. The common components are the 
uncertainty in the measured integrated luminosity, the 
uncertainty in tracking and the uncertainty in particle 
identification. Since the absolute branching fraction scale 
for D + — * if ~7r + 7r + and D° — > K~Tr + ir + ir~ depend on 
the branching fraction scale for D° — * K~ir + , the to- 
tal percentage uncertainty for the two channel branching 
fractions (6.6% and 4.1%) are split into a common com- 
ponent that matches the percentage uncertainty for the 
D° — > if ~7r + branching fraction (2.3%) and indepen- 
dent components (6.2% and 3.4%). All other systematic 
uncertainties are treated as independent and added in 
quadrature. The common uncertainties are added lin- 
early. 

As a comparison, Table IV lists the observed cross sec- 
tions for D°D° and D + D~ production at the cm. ener- 
gies of 3.773 GeV and 3.768 GeV, which were measured 
by this experiment and MARK-III [Hi ■ 



VI. RADIATIVE CORRECTIONS 

In any e + e~ colliding beam experiment, the electron 
(positron) always radiates at the interaction point be- 
cause of the potential of the positron (electron). Since 
this radiation (Bremsstrahlung) carries energy away, 
the actual center-of-mass energy for the e + e~ annihila- 
tion is reduced by Bremsstrahlung to y/s(l — x), where 
£i?bcam is the total energy of the emitted photons. The 
Bremsstrahlung is principally responsible for the distor- 
tions to the tree level resonance line shape, while the self 
energy of the electron and positron and the vertex cor- 
rections to the initial state affect the overall factors to 
change the scale of the cross section. All of these correc- 
tions are called Initial State Radiation (TSR.~) corrections. 
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The tree level cross section for DD production at the 
energy of 3.773 GeV can be obtained by correcting the 
observed cross section for the effects of the ISR and vac- 
uum polarization. 

The observed cross section, cr obs , at the nominal energy 
y/s can be written as a convolution of the Born cross 
section o B (s(l — x)) and a sampling function f(x, s), 



5? + 



(14) 



(15) 



a ohs (s) = / dx-f(x,s)a B {s(l-x))(l + 6 VP (s(l-x))). 
Jo 

The vacuum polarization correction (1 + 5yp) includes 
both leptonic and hadronic terms. It varies from charm 
threshold to 4.14 GeV by less than ±2% Q. In this data 
analysis, we treat it as a constant of 



(1 + Svp) = 1.047 ± 0.024. 



(8) 



Since we are interested in the ^>(3770) resonance in this 
analysis, we take the a B to be the bare Breit-Wigner 
cross section 



\E) 



i27rr° e r tot (£) 



(9) 



where T° ee = T ee / (1 + <5 vp ), M and T ee are the mass and 
leptonic width of the ^(3770) resonance respectively; E is 
the center-of-mass energy; Y to t(E) is chosen to be energy 
dependent and normalized to the total width r to t at the 
peak of the resonance Q • The r to t (E) is defined as 



PnO 



r to t(£) = r 



l + (rp D o) 2 l+(rp D ±) 2 
P nO , P D± 



(10) 



l + (rp° D o) 2 ^ l+(rp° D ±) 2 



where p° D is the momentum of the D mesons produced 
at the peak of "0(3770), p D is the momentum of the D 
mesons produced at the cm. energy yfs, Tq is the width 
of the -0(3770) at the peak, and r is the interaction radius 
which was set to be 0.5 fm. In the calculation of the Born 
order cross section, the 0(3770) resonance parameters 
M = 3769.9 ± 2.5 MeV; T = 23.6 ± 2.7 MeV and T ee = 
0.26 ± 0.04 keV @ were used. 

In the structure function approach introduced by Ku- 
raev and Fadin [12] 13], the sampling function can be 
written as 



f(x,s) = I3x p - 1 5 v+s + S H , (11) 
where [3 is the electron equivalent radiator thickness, 



0= -(^4 



1 



mf. 



(12) 



3 a /tt 2 1 

= 1+ i^ y-2 



fl 2 



24 V 3 



] m4 



2tt z - 



37 



4 



6? = 



8' 



4(2 - x) In - - 

x 



1 l + 3(l-a;) 5 



ln(l-x) -6 + 
(16) 

In the above formula, m e is the electron mass and a is the 
fine structure constant. The 0(3770) width (~ 24 MeV) 
is much large than the energy spread (~ 1.37 MeV) of 
BEPC. So the effect of the beam energy spread on the 
cross section could be ignored. The 0(3770) is gener- 
ally assumed to decay exclusively into DD. Taking these 
considerations, the observed cross section of Equation (7) 
should be replaced by 



obs 



(s) = (1 + S VP ) 



l-4MJ/s 



dx ■ f{x,s)a B (s(l- x)) 
(17) 



in calculation of the radiative corrections. The correction 
factor for the radiative effects is given by 



T obs 



9 



(18) 



Figure 2 shows the factor of the radiative corrections as 
a function of the nominal center-of-mass energy. At the 
center-of-mass energy y/s = 3.773 GeV, the factor is 



g = 0.779 ±0.031, 



(19) 



where the error is the uncertainty arising from the errors 
on the ^(3770) resonance parameters. The uncertainty 
is mainly due to the error of the mass of the resonance. 



VII. CROSS SECTION FOR DD PRODUCTION 

The tree level cross section for DD production is ob- 
tained by dividing the observed cross section by the fac- 
tor of the radiative corrections. At y/s = 3.773 GeV, the 
charged, neutral and total tree level D pair production 
cross sections are 



a D o BO = (4.60 ± 0.12 ± 0.45) nb, 



a D+D - = (3.29 ±0.10 ±0.37) nb, 



and 



'DD 



(7.88 ±0.15 ±0.74) nb, 



(20) 



(21) 



(22) 



where the first error is statistical and the second system- 
atic which include the uncertainty in the factor of the 
'radiative corrections. These results are compared to the 
connled-cha.nnel model nrediction in Table V. 
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3.76 3.78 3.8 3.82 3.84 3.86 3.88 3.9 



Nominal E r 



[GeV] 



FIG. 2: The factor of radiative corrections as a function of 
the nominal center-of-mass energy. 



TABLE V: Comparison of tree level cross section measure- 
ments with prediction of the coupled-channel model at y/s = 
3.773 GeV. 





Experiment 


coupled-channel Model 


a D o£,o 


4.60 ±0.12 ±0.45 nb 


1.80 nb 


°D+D- 


3.29 ±0.10 ±0.37 nb 


1.28 nb 




7.88 ±0.15 ±0.74 nb 


3.08 nb 



BES-II experiment measured R u ds 01 > which is the ratio 
of the tree level light hadron (containing the u,d and s 
light quarks) cross section over that for production 
in the energy region from 2.0 to 3.0 GeV. Theoretical 
expectation is that R uc i s is approximately independent 
of center-of-mass energy in this region [3j. Fitting to 
the R u ds values at 9 energy points in the energy region, 
we obtain R uds = 2.26 ± 0.14. Assuming that 0(3770) 
decays exclusively into DD, the value of R is evaluated 
using R — i?£>/2 ± R u ds, which gives 



R= 3.55 ±0.03 ±0.18. 



IX. SUMMARY 



(25) 



In summary, using the 17.3 pb _1 of data collected with 
the BES-II detector at BEPC at center-of-mass energy 
\/i = 3.773 GeV, the observed cross sections for D°D°, 
D + D~ and DD production have been measured. Those 



are a 



D°D" 



(3.58 ±0.09 ±0.31) 

T obs 



nb, 



T obs 
T D+D- 



(2.56 ± 



0.08 ± 0.26) nb and a°° s D = (6.14 ± 0.12 ± 0.50) nb. The 
tree level cross sections for the D°D°, D + D~ and DD 
production are determined to be (T d oqo = (4.60 ± 0.12 ± 
0.45) nb, (J D+D - = (3.29 ± 0.10 ± 0.37) nb and <j dD = 
(7.88±0.15±0.74) nb, which are about a factor 2.5 times 
larger than that predicted by the coupled-channel model. 
Using the measured R u ds i n the energy region from 2.0 
to 3.0 GeV from BES-II experiment and assuming that 
-0(3770) only decays to DD, the total tree level cross 
section for inclusive hadronic event production at 3.773 
GeV is obtained to be R = 3.55 ± 0.03 ± 0.18. 



VIII. MEASUREMENT OF R D AND R 
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